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Evidence of first-order transition between vortex glass and Bragg glass phases 
in high-T c superconductors with point pins: Monte Carlo simulations 

Yoshihiko NonomuraQ and Xiao Hu 
National Research Institute for Metals, Tsukuba, Ibaraki 305-0047, Japan 
(Received 14 February 2000) 

Phase transition between the vortex glass and the Bragg glass phases in high-Tc superconductors in 
B || cis studied by Monte Carlo simulations in the presence of point pins. A finite latent heat and a 
5-function peak of the specific heat are observed, which clearly indicates that this is a thermodynamic 
first-order phase transition. Values of the entropy jump and the Lindemann number are consistent 
with those of melting transitions. A large jump of the inter-layer phase difference is consistent with 
the recent Josephson plasma resonance experiment of P^S^CaC^Os+j, by Gaifullin et al. 
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Vortex states in high-T c superconductors have been 
intensively studied experimentally and theoretically |Q. 
Because of large fluctuations owing to high transition 
temperature and strong anisotropy, the flux-line lattice 
(FLL) melts at much lower temperatures than those pre- 
dicted by Abrikosov's mean-field theory. The FLL melt- 
ing is a thermodynamic first-order phase transition. In 
pure systems, the melting line stretches up to a high mag- 
netic field as large as H C2 . However, all experiments show 
that first-order melting lines terminate at much lower 
magnetic fields ||||. Complicated phase diagrams are 
obtained experimentally in Bi2Sr2CaCu20s+ a (BSCCO) 
|| and YBa 2 Cu 3 7 _ 5 (YBCO) §, and it is believed that 
effects of impurities are essential in real materials. For 
example, vacancies of oxygen atoms, which play the role 
of point pins to flux lines, cannot be excluded completely 
even in crystals of highest quality. 
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FIG. 1. Schematic vortex phase diagram of high-Tc su- 
perconductors with point pins. 



Fisher et al. || studied the Ginzburg-Landau model 
in a random potential, and proposed the so-called vortex 
glass (VG) phase before experimental studies. Giamarchi 
and Doussal pointed out that the Bragg glass (BG) 
phase can exist in weak fields. In this phase, correlation 
functions decay in power laws |^,^|, and the structure 
factor shows a triangular Bragg pattern. Accordingly, 
a phase transition between these two glass phases may 
be observed when the magnetic field is sweeped across 
the phase boundary (see Fig. |l|). The existence of the 
VG-BG phase transition in vortex systems with point 
pins shows a sharp contrast to the phase diagram of pure 
systems. Presuming a first-order phase transition, the 
VG-BG phase boundary was evaluated phenomenologi- 
cally (h]. ll|]. However, physical properties around this 
phase boundary have not been clarified in experiments 
and numerical calculations until recently. The shape of 
the phase boundary seems to depend on observed quan- 
tities, and some experiments even suggest a crossover 
rather than a phase transition. Although some simula- 
tions [p^-|l5t gave phase diagrams similar to that of Gia- 
marchi and Doussal, numerical accuracy of these studies 
was not good enough to distinguish phase transitions and 
crossovers. The stability of the VG phase was studied by 
Kawamura Jpc| including the screening effect. 

Quite recently, Gaifullin et al. observed [[l7| a large 
jump of the inter-layer phase difference on the VG-BG 
phase boundary of BSCCO by the Josephson plasma res- 
onance experiment. They claimed that their observation 
is the evidence of a first-order phase transition. In the 
present Letter, we show more direct evidence of the first- 
order phase transition on the VG-BG boundary by large- 
scale Monte Carlo simulations. That is, a finite latent 
heat and a (^-function peak of the specific heat are ob- 
served. Sharp jumps of the inter-layer phase difference 
and the averaged fluctuations of flux lines are also ob- 
tained. 

In order to clarify vortex states and phase transitions of 
high-T'c superconductors in the presence of point pins, we 
start from the three-dimensional anisotropic, frustrated 
XY model on a simple cubic lattice [Q[L9|. Effects of 
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point pins are introduced into the model by randomly- 
distributed weakly-coupled plaquettes in the ab plane. 
Since a vortex sitting on a plaquette costs an energy pro- 
portional to the couplings surrounding it, flux lines tend 
to penetrate plaquettes with weaker couplings in order 
to reduce such loss of energies. The Hamiltonian of our 
model is given by 

H = - ^2 Jij cos (<fo - fa - Aij) 

i.j£ab plane 

-p2 COS ((j) m - <j) n ) , (1) 

m,n || c axis 

Ay ' = Sf A(2) ' dr(2) ' (2) 

with the periodic boundary condition along all the direc- 
tions. Couplings in the ab plane are given by Jy = bJ 
(0 < b < 1) on the weakly-coupled plaquettes, and 
Jij = J otherwise. The density and the strength of point 
pins are controlled by the probability of weakly-coupled 
plaquettes, p, and the parameter b, respectively (see Fig. 

. The pinning energy is of order of (1 — b) J. A uni- 
form magnetic field is applied along the c axis, and its 
strength is proportional to the averaged number of flux 
lines per plaquette, /. Here we concentrate on the model 
with L x — L y — 50 and L c = 40. This system size is 
large enough to describe the melting transition in the 
pure system (6 = 1) p9| . 
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FIG. 2. Schematic description of point pins in the present 
model. 

In our model, we have four adjustable parameters: the 
anisotropy constant T, the density of flux lines /, the 
density of point pins p, and the strength of point pinning 
b. In order to investigate the VG-BG transition, we vary 
6, while fix the temperature at T = 0.06 J//cb and other 
parameters at T = 20, / = 1/25 and p = 0.003. In other 
words, material parameters of the bulk system and the 
number and positions of point pins are not changed dur- 
ing the simulations. As will be shown later, this temper- 
ature is low enough for the study of the VG-BG phase 



boundary. Typical Monte Carlo steps (MCS) with the 
Metropolis algorithm are 3~4x 10 7 MCS for equilibra- 
tion, and 0.5 ~lx 10 7 MCS for measurement. 
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FIG. 3. Internal energy e (squares) and specific heat C 
(triangles) per flux line per ab plane versus strength of point 
pinning b. Straight lines are drawn as guides for eyes. 

First, the b dependence of the internal energy e and 
the specific heat C per flux line per ab plane is displayed 
in Fig. ||. Clearly, the internal energy shows a sharp 
jump at the transition point b* = 0.895 ± 0.005 with a 
latent heat per flux line per ab plane Q w 2.3 x 10~ 2 J, 
and a 5-function peak of the specific heat also occurs at 
the same parameter. These two facts indicate that the 
VG-BG transition is a thermodynamic first-order phase 
transition. From this latent heat, the entropy jump at b* 
is estimated as 

AS — Q/T w 0.38fc B , (3) 

which is comparable to the experimental value in the 
melting transition of YBCO, AS « 0.5/fc B |o|- 

Second, the b dependence of the inter-layer phase dif- 
ference, (cos(<^„ — <j) n +i )), is plotted in Fig. [|. This 
quantity is related to the Josephson energy per phase 
variable ej and the anisotropy constant T as 

(cos(0„ — (pn+i)) = -e,jr 2 /J, (4) 

and a small change of ej is magnified in this quantity 
in extremely anisotropic systems. This quantity also 
jumps sharply at b* — 0.895 ± 0.005, and the value of 
the jump at b* , Apd ~ 0.12, is as large as the experi- 
mental value, Apd ~ 0.2 ]17[] . Moreover, the ratio of the 
jump of the Josephson energy to the latent heat is given 
by Aej/(Q/) w 0.34, which means that the latent heat 
is equally distributed to all the directions in the VG-BG 
phase transition in extremely anisotropic systems. 
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FIG. 4. Inter-layer phase difference (cos(<j!>„ — 4>n+i)) ver- 
sus strength of point pinning b. Straight lines are drawn as 
guides for eyes. 

Third, the Lindemann number is evaluated directly 
pl| . The deviation u of a flux line is measured in each 
ab plane from the projection of the mass center of the 
flux line, and averaged over all the flux lines and the ab 
planes. Then, the Lindemann number cl is given by 



cl = lim {u 



2x1/2 



/«0 , 



(5) 



with the lattice constant of the triangular FLL, do = 
(2/ v / 3) 1/2 // 1/2 . The b dependence of {u 2 ) 1/2 /a Q is 
shown in Fig. H and we have cl ~ 0.28. This value 
is almost equal to the one obtained in the FLL melting 
of pure systems plj . 
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FIG. 5. Averaged fluctuations of flux lines (u 2 ) 1 ' 2 /ao ver- 
sus strength of point pinning b. 



Finally, we go into some details of the present simu- 
lations. The system with b = 0.90 is calculated at first. 
Simulations are started from a random configuration at 
a very high temperature, and the system is gradually 
cooled down to T = 0.06J//cb- During the cooling pro- 
cess, the first-order melting transition characterized by a 
discontinuous appearance of the helicity modulus along 
the c axis, T c |jl9), takes place at T m 0.079 J/fcs, which 
corresponds to the vortex liquid (VL)-BG phase transi- 
tion. Then, the strength of point pinning b is varied. 
Since the quantity T c is proportional to the superfluid 
density, the region with finite T c is superconducting. 
This quantity is nonvanishing for all the values of b shown 
in Figs. ||-^| at T = 0.06J/fce, an d therefore the phase 
transition investigated in the present Letter is not the 
VL-BG one, but the VG-BG one. Equilibration in sys- 
tems with point pins is much slower than that in pure 
systems, and only one sample can be taken for calcula- 
tions at present. Nevertheless, the results obtained in the 
present Letter are quite clear-cut and consistent with ex- 
periments. Thus, the small number of random sampling 
does not seem serious. Since positions of point pins are 
independent in each ab plane, the number of ab planes, 
L c = 40, would be large enough for averaging effects of 
point pins. 

Although we have concentrated on the VG-BG tran- 
sition for a single density of point pins p in the present 
Letter, we have also investigated the VL-BG and VL-VG 
transitions for various p, and obtained the overall phase 
diagram in the p-T plane. The structure of the p-T 
phase diagram is similar to that of the B-T phase dia- 
gram. Experimentally, the increase of p corresponds to 
the repeated irradiation of electrons or protons, and our 
p-T phase diagram is consistent with recent experiments 
[22,E3[. Details of this study will be reported elsewhere 

A 

In conclusion, the first thermodynamic evidence of the 
first-order transition between the vortex glass (VG) and 
the Bragg glass (BG) phases has been obtained in high- 
T c superconductors in the presence of point pins. A finite 
latent heat and a <5-function peak of the specific heat are 
observed by large-scale Monte Carlo simulations of the 
three-dimensional anisotropic, frustrated XY model with 
randomly-distributed weakly-coupled plaquettes. The 
entropy jump derived from the latent heat is nearly equal 
to those in the melting transition of YBCO. The Linde- 
mann number evaluated from fluctuations of flux lines, 
cl ~ 0.28, is reasonable for the first-order phase transi- 
tion. The inter-layer phase difference also shows a sharp 
jump on the VG-BG phase boundary. This property is 
consistent with the Josephson plasma resonance experi- 
ment of BSCCO by Gaifullin et al. 
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